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Abstract 

In order to check on a recent suggestion that local scale invariance [M.Henkel et al. 
Phys. Rev. Lett. 87, 265701 (2001)] might hold when the dynamics is of Gaussian nature, we 
have carried out the measurement of the response function in the kinetic Ising model with Glauber 
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t^ \ dynamics quenched to Tq in d = 4, where Gaussian behavior is expected to apply, and in the two 

other cases of the d = 2 model quenched to Tq and to below Tq, where instead deviations from 

o ; 

Gaussian behavior are expected to appear. We find that in the d = 4 case there is an excellent 

(N '. 

agreement between the numerical data, the local scale invariance prediction and the analytical 

fN| ' Gaussian approximation. No logarithmic corrections are numerically detected. Conversely, in the 

en ' 

d = 2 cases, both in the quench to Tq and to below Tq, sizable deviations of the local scale invari- 
O , 

ance behavior from the numerical data are observed. These results do support the idea that local 
O 
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scale invariance might miss to capture the non Gaussian features of the dynamics. The considerable 
precision needed for the comparison has been achieved through the use of a fast new algorithm for 
the measurement of the response function without applying the external field. From these high 



T3 

a 

q I quality data we obtain a = 0.27 ± 0.002 for the scaling exponent of the response function in the 



d = 2 Ising model quenched to below Tq, in agreement with previous results. 
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I. INTRODUCTION 

The non equilibrium dynamics of aging and slowly evolving systems is a topic of current 
and wide interest [l|. Much work has been devoted to the understanding of two times 
quantities such as the auto-correlation function C(t, s) = (<j)(x,t)<j)(x,s)) and the auto- 
response function R(t,s) = 5(<f)(x,t)) /5h(x, s), where <fr(x,t) is the order parameter at the 
space-time point (x, t), h(x, t) is the conjugate external field and the averages are taken over 
the thermal noise and the initial condition with t > s. The interest in the relation between 
these two quantities dates back to the solution by Cugliandolo and KurchanJ2| of the p-spin 
spherical model, where they introduced the fluctuation-dissipation relation as a measure of 
the distance from equilibrium. Furthermore, this relation can encode important information 
on the structure of the equilibrium state J2| • 

One among the simplest examples of systems exhibiting aging and slow dynamics is a 
ferromagnetic model evolving with a dissipative dynamics after a quench from an infinite 
temperature to a final temperature T smaller than or equal to the critical temperature Ta- 
in both cases, the slow relaxation entails the separation of the time scales. That is, when s 
becomes large enough, the range of r = t — s can be devided into the short rCs and the long 
t 3> s time separation, with quite different behaviors in the two regimes. The first one is the 
quasi-equilibrium or stationary regime, where the two time quantities are time translation 
invariant (TTI) and exhibit the same behavior as if equilibrium at the final temperature of 
the quench had been reached. The second one is a genuine off equilibrium regime, where 
aging becomes manifest. A crucial point is how these two behaviors are matched. The generic 
pattern for phase-ordering systems is that the matching is multiplicative in the quences to 
Tc and additive in the quenches to below Tc- The first one is well documented by analytical 
calculations. For the second one, although the analytical evidence is less abundant, the 
additivity is required on general grounds by the weak ergodicity breaking scenario |l|. 

To be more specific, in the case of the quench to Tc, using the methods of the field 



theoretical renormalization group (RG), the evolution equations for C(t,s) aiKLR 
obtained by means of a series expansion around the Gaussian fixed point 



1, s) are 
The 



solution of these equations gives for R(t, s) the scaling form 

R(t,s) = s-^F R (t/s) (1) 

with the additional requirement that the scaling function must be of the form 

F R (x) = A R (x - l)^ 1+ ^x e f R {x) (2) 

where a = {d — 2 + r])/z, d is the space dimensionality, r] and z are the usual static and 
dynamic critical exponents, 9 is the initial slip exponent and lmx^oo f R (x) = 1 J4J, |5j. A 
similar result is obtained for the correlation function |4j, |5|, |yj. The multiplicative structure 
becomes evident rewriting Eq. (JTJ) as 

R(t lS ) = A R (t-s)-^g R (x) (3) 

where g R (x) = x e f R (x). 

In the case of the quench to T < Tc, the above form is replaced by the additive structure 

R{t,s) = R st {t-s) + R ag {t,s) (4) 

where the first is the stationary contribution and the second one is the aging contribution 
which obeys a scaling form of the type ((TJ) 

Ra 9 (t,s) = s^ 1+a ^h R (x) (5) 

without the restriction (J2J) on the form of h R (x). 

In the latter case, all theoretical efforts have been direct toward the determination of 
the exponent a and the scaling function h R (x). Keeping into account that the evolution is 
controlled by the T = fixed point, which in no case is Gaussian, the perturbative RG cannot 
be used and resorting to uncontrolled approximations is unavoidable. Among these, one of 
the most successful is the Gaussian auxiliary field (GAF) approximation. The method was 
originally introduced by Ohta, Jasnow and Kawasaki jj| in the study of the scaling behavior 
of the structure factor and has been subsequently applied also to the study of the response 
function [alsl]- Recently, new results for R(t,s) have been obtained by Mazenko [l0| using 
a perturbative expansion which improves on the GAF approximation. Next to approximate 
methods, there exist exact analytical results for two solvable models: the one dimensional 



Ising model 



Q 



12j and the O(N) model in the large A^ limit for arbitrary dimensionality 



[13[. Both solutions give for R ag {t,s) the scaling form (jSJ), with a = for the d — 1 Ising 
model and a = (d — 2)/2 for the large iV model. 

Therefore, in the context of the quenches to below Tc, where a controlled theory is not 
available and numerical simulations are very time demanding, it is of much interest the 
CO**™ put forward by He„ M et a, fi Q that the _ taction transform, 
covariantly under the group of local scale transformations, both in the quenches to and to 
below Tc- The hypothesis of local scale invariance (LSI), then, implies that the multiplicative 
structure for R(t, s), as obtained from RG arguments at T c , applies also in the quenches to 
below Tq- That is, from LSI follows that R(t, s) obeys Eq. (J3J), both at and below Tc, with 
the additional prediction that 

f R (x) = 1 (6) 



holds not just asymptotically, but for all values of x, while the amplitude Ar and the 
exponents a and 9 remain unspecified. Hence, whith the LSI hypothesis, the difference 
between the quenches to Tc and to below Tc would be left only in the values of the exponents 
a and 9. This is actually verified by the exact solution of the spherical model [Lj, llfi| . 
Conversely, from the GAF approximation and from the exact solution of the d = 1 Ising 
model follows 17[ 

f R {x) = [(z-l)/z] 1 '' (7) 

which differs significantly from the above LSI prediction ©. 

In the case of the quench to T = Tc, the validity of LSI has been tested by Calabrese 
and Gambassi |l8| by means of the e expansion. Their field theoretical computation shows 
that LSI holds up to the first order in e = 4 — d, but deviations of order e 2 are present. 
Motivated by this result, Pleimling and Gambassi (PG) in a recent paper [l9| have carried 
out a careful numerical check of both LSI based and field theoretical calculations in the 
Ising model quenched to T c , in d = 2 and d = 3. In particular, they have computed the 
integrated global response to a uniform external field, finding i) a discrepancy between the 
LSI behavior and the data, ii) that the discrepancy is more severe in d = 2 than in d = 3 
and iii) that the e 2 correction does not eliminate the discrepancy, but improves on the LSI 
prediction. In this connection, Calabrese and Gambassi [y] first and then PG made the 
remark that the LSI prediction coincides with the Gaussian approximation, thus accounting 
for the agreement between LSI and the solution of the spherical model. 



Following through this suggestion, one could anticipate that the discrepancy between 
the LSI behavior and the results of simulations should disappear in the quench to Tc with 
d — 4, while it ought to get even worse in the quench to below Tc, independently of the 
dimensionality. Furthermore, in the latter case the failure of LSI is expected to be not just 
in the quantitative accuracy of the approximation, but also of a structural character since 
the multiplicative form of R(t, s) is incompatible with the weak ergodicity breaking scenario. 

In order to investigate these ideas, one can take advantage of the efficient numerical 

nnfl 

tools made available by a new generation of algorithms |2(J, |21|, |22J. These algorithms are 



based on the relation between R(t, s) and unperturbed quantities which, by speeding up the 
or the measurement of R(t,s). In this paper, exploiting the algorithm 



simulation, allow : 

introduced by us 22j, we extend the investigation of the Ising model carried out by PG 
to the two cases of the quenches to Tc with d = 4 and to below Tc with d = 2. Rather 
than computing the integrated response function for a global quantity, as PG have done, we 
access directly the local response function R(t, s), thus making the comparison between the 
numerical data and Eq.fjHJ). 

In the d = 4 Ising model quenched to T = Tc, after addressing the question of the univer- 
sality of the exponent 9 23| and of the ratio TcAr/Ac between the amplitudes of response 



and correlation function 



i III llii |24j , we find an excellent quantitative agreement between 
the numerical data and the analytical results from the Gaussian model. In particular, we 
find that both for R(t, s) and C(t, s) not only the scaling exponents, but also the scaling 
functions and the ratio TcAr/Ac are well accounted for in the Gaussian approximation. We 
find that Eq. @ holds and we conclude that LSI correctly describes the critical quench of 
the d — 4 Ising model. Conversely, in the quench of the d — 2 Ising model to T = Tc and 
to T < Tc, important deviations from LSI are observed. These findings do bring support 
to the idea that the LSI principle is some sort of zero order theory of Gaussian nature and 
contradict previous statements jlflL l25| that no deviations from LSI predictions are observed 
in the measurements of local quantites. 

The paper is organized as follows. In sec]TI]we shortly review existing results for Cit, s) 
and R(t, s). In particular in sec JII Al and in sec lIIBI we give the results from RG arguments 
and from the Gaussian model, respectively, while in sec III CI we present a phenomenological 
picture for the quench to T < Tc- In sec JIHI we outline the algorithm used in the simulations 
and in sec II VI we present and discuss the numerical results. Concluding remarks are made 



in the last section. 



II. EXISTING RESULTS 



We consider a system with a non conserved scalar order parameter 0(af, t) (model A in 
the classification of Hohenberg and Halperin 26[) evolving with the Langevin equation 

d(j)(x,t) 8H[<j>] 



dt 8<f>{x,t) ' V ^ 

where i](x,t) is a Gaussian white noise with expectations 

( v (x : t)) = (n(x, t)r](x', t')} = 2T5(x - x')5(t - t') 

and H[(j)\ is of the Ginzburg-Landau- Wilson form 

1 



(8) 



(9) 



H[<t>] 



dx 



_(V0) 2 + Ir0 2 + i^ 4 



(10) 



with r < and g > 0. The system is prepared in an uncorrelated Gaussian initial state with 
expectations 

(0(f , 0)) = (cf>(x, O)0(f , 0)) = r G l 5{x - x'). (11) 



A. Quench to Tq: RG results 

In the case of the quench to Tc one can show, by means of standard RG methods |4|, 
that Tq" 1 is an irrelevant variable. Thus, putting Tq" 1 = 0, one obtains the leading scaling 
behavior which is given by Eqs. (|1I2|) for R(t,s), whereas for the correlation function one 
has 

C(t,s) = s- b F c (t/s), (12) 



with the scaling function 



F c (x) = A c {x - l)- b x y - L f c (x) 



and 



b = a 



d-2 + f] 



(13) 



(14) 



As for Jr(x), the RG method allows to fix only the large x behavior linx^oo fc(x) = 1. 



From Eq.(|14|) one has that a and b are related to the critical exponents rj and z. Therefore, 

2611. a and b take the same value 



according to the classification of Hohenberg and Halperin 

for systems belonging to the same class of universality. The problem of the univerality of 9 

has been addressed in a series of papers J23l. l24j| . Furthermore, Godreche and Luck Il6j have 



proposed that also the ratio TqA^/Ac is a universal quantity More precisely, considering 



the limit fluctuation dissipation ratio [2j X^ defined by 



T-l 

ioQ 



X^ = ton ton J^|4 (15) 

s^oct->oo d s C(t,S) 

and using Eqs. ()ll2ll2|) . in the quench to the critical point one has 

Universality of 9 and TcAr/Ac implies universality of X^. Indeed, we will see that numer- 
ical results for 9 and X^, in the d = 4 Ising model, give the same values as in the Gaussian 
model. 



B. Quench to Tq- the Gaussian model 

The critical Gaussian model is obtained putting r = and g = in the Hamiltonian (jl(Jj) . 
Then, the equation of motion (|HJ) can be solved in Fourier space yielding 



T, 



k 



C(k, t, s) = (0(fc, t)4>(-k, s)) = 4 e~ k2 ^ - e~ k2 ^ + (17) 



e 



k\t+s) 



T(> 



#,M) = » = e^M (18) 

Sh(—k, s) 

whith t > s. The auto-correlation function and the auto-response function are obtained 
integrating over k the above equations. In order to regularize the equal time behavior of 
C(t, s) and R(t, s) one must introduce an high momentum cut-off and, for simplicity, we 
choose a smooth cut-off implemented by the multiplicative factor e~ k / A in Eqs. (|17ll8J) . 
Neglecting the last term in Eq.(|17|). in order to keep only the leading scaling behavior, one 
gets |22| 

C(t, s) = {d _ 2) ^ n)d/2 [{t-* + to) 1 -** -{t + s + to) 1 -* 2 ] , (19) 

R(t,8) = j^(t- 8 + *)-** (20) 
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where to = 1/A 2 . Notice that the specific choice of the cut-off affects the behavior of R(t, s) 
and C(t, s) only on the time scale t — s ~ to- Taking t — s ^> to, the above results are in the 
scaling form of Eqs. (|13I2|) with fn(x) = 1, as required by LSI, and with a = b = d/2 — 1, 
# = 0, fc{ x ) — x ~ x(x — l) a (x + l) _a . In particular, in d = 4 one has 

C(t, s) = A c s~\x - 1 + t / S y\x + 1 + to/s)- 1 (21) 

i2(f,s) =A iJ (t-s + t )- 2 , (22) 

with A c = 2T c /(4vr) 2 and 2T C A R = A c . 

C. Quench to below Tc- phenomenological picture 

In the case of the quench to below Tc, the system evolution is characterized by the 
formation and subsequent growth of compact ordered domains whose typical size increases 
with the power law 

L(t) ~ t 1/z . (23) 

The evolution via domain coarsening suggests [l|, for large s, the additive form of the 
correlation function 

C(t : s) = C st (t-s) + C ag (t,s) (24) 

where C st (t — s) represents the correlation function of the equilibrium fluctuations within an 



infinite domain and C ag (t, s 
as well as numerical results 



13 



16] 



is the domain walls contribution. Analytical solutions 
28| confirm this structure, with C ag (t, s) obeying a scaling form 
as in Eq.(|12j) and with b = 0. As stated in the Introduction, the similar structure (@J) holds 
also for the response function with R ag {t, s) in the scaling form (jHJ) and R s t{t — s) related to 
C s t{t — s) by the fluctuation dissipation theorem, TR st (t — s) = dC st (t — s)/ds. 

Numerical simulations [9|, |l7( for the zero field cooled magnetization x(t, t w ) = 
J. R(t, s)ds are consistent with the additive structure (JU) and with a scaling function in 
Eq. (J5J) of the form 

Mx) = AR ( x -i + l/sy-v*» (25) 

where Ar, a, (3, to are phenomenological parameters, while z is the dynamical exponent en- 
tering Eq. (|23|). Here, t is a microscopic time which is negligible except when x — > 1. Recent 
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results |28( from the direct measurement of R(t, s) do support the above form of Hr{x) and 
give a quantitative estimate of a, Ar and j3. The physical meaning of Eq. (J25J1 becomes 
clear for short time separation t — s <C s. In this case Eqs. (|5I25|) can be rewritten as 

Rag(t, S) = pl(s)R sing (t - S) (26) 

where 

Rsing(t ~S)= A R (t - S + to)- 1 - * 1 '' (27) 

and pi(s) oc L~ l (s) is the interface density at time s. Therefore, R S i ng (t — s) represents 



the response of a single interface and Eq. (jzbj) simply means that the aging contribution in 
the response is produced by the interfacial degrees of freedom. For larger time separation 
t — s > s, interfaces interact with each other and the interaction generates the term x 13 
in Eq. (|25j). The form ()25|) of Hr{x) is corroborated by the exact analytical result for the 
d = 1 Ising model with non conserved order parameter jlli Il2j |. by the numerical results 
of the d — 1 Ising model with conserved order parameter [22| and by the analytical results 
obtained with the GAF approximation [8|, |9|. In the sec lIVDl we will present a direct 
comparison between Eq. (|25|) and the prediction from LSI. 

III. THE ALGORITHM 

We consider a system of N spins on a lattice with the Ising Hamiltonian 

H = -jY] OiO d (28) 



E 



where the sum runs over the nearest neighbours pairs < ij > and J > 0. The time evolution 
is then obtained through single spin flip dynamics with Glauber transition rates 



Wi([a] -»■ [a']) = - 



h w 



(29) 



where [a] and [a'} are spin configurations differing only for the value of the spin in the i- 
th site, hY = J ^2 k <7k is the Weiss field, J is the ferromagnetic coupling and the sum is 
restricted to the nearest neighbors of the z'-th site. C(t, s) and R(t, s) are given by 

CM = ^E^(*)^( S )> (30) 



and 



R(t, s) 



lim 

As^O AsN 






d{a l {t)) 



dh. 



(31) 



h=0 



where &i(t) is the spin in the i-th site at time t and hi is an external field acting on the 
i-th site during the time interval Is, s + As]. In the computation of R(t,s), we use our 

n nn 

own algorithm [22j], which offers higher efficiency with respect to other methods [20, 12 lj 
allowing to compute the response function without imposing the external field. Carring out 
the derivative in Eq. (|31j) we find 



TR{t, s) 



- lim 

2 As^O 



C(t, s + As)- C(t, s) 
As 



(ai(t - As)Bi(s)) 



where Bi enters the evolution of the magnetization according to 

d{<n(t)) 



22J 



dt 



mt)). 



(32) 



(33) 



The above result is quite general and is independent of the details of the Hamiltonian and 
of the transition rates. Furthermore, it can be easily generalized to the case of vector order 
parameter |2j| . In the case of single spin flip dynamics, one has 



Bi(t) = 2a i {t)w i {[cr] -> [<x { ]) 



(34) 



with Wj([cr] —* [cTj]) given in Eq.( 

In order to improve the signal to noise ratio, we compute the quantity 



fi,(t,[s+l,s]) 



s+l 



R(t,t')dt' 



(35) 



which is the response to a perturbation acting in the time window [s, s + l]. Here and in the 
following we express time in units of a Monte Carlo step. Replacing the integral in Eq. (J35|) 
by a discrete summation on the microscopic time scale e = 1/N, from Eq. (|3*2*|) we obtain 



T/i(£,[ S + M]) = - 



1 N 
C(t, s + l) - C(t, s)-jjJ2 ( a ^ - l l N ) B ^ s + ( k - iVAO) 



i,fe=i 



• (36) 



Because of the scaling form (JIJ) for R(t, s), it is easy to show [28| that R(t, s) coincides with 
fi(t, [s + l,s]) up to corrections of order 1/s which, in the considered range of times, can 
always be neglected. Therefore, in the simulations we identify R(t, s) with //(£, [s + l, s\) and 
the numerical results for R(t, s) are obtained from Eq. (|36|) . In all cases we take a completely 
disordered initial state which, in principle, produces a correction to scaling. However, this 
is not detectable in the time region explored in the simulation. 
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IV. NUMERICAL RESULTS 



d = 4, T = T c 



We have considered a system of iV = 60 4 Ising spins on a four dimensional hypercubic 
lattice quenched to the critical temperature Tc — 6. 68 J |30j. The response and correlation 
functions are then computed for four different values of s = 25, 50, 75, 100. In all Figures 
the error bars are smaller than the symbols. 



In order to compare with the results of the Gaussian model given in III Bl we observe 
that R(t, s) in Eq. (J22)) depends only on the time difference t — s. This result is reproduced 
in the numerical simulations. Indeed, plotting the curves for different s as a function of 
t — s (see Fig C] we find the collapse on a master curve that is well described by Eq. (|22|). 
with TqAr = 0.35±0.02 and to — 0.1. There is a very small difference only for short time 
separations t — s ~ 1, which can be attributed to the specific choice of a smooth cut-off used 
in the integration over k of Eq. IjIB]) . 

In Fig|21 we compare the numerical data for C(t, s) with the analytical expression of 
Eq. (J21JI . The plot shows that the quantity sC(t, s) depends only on the ratio t/s in complete 
agreement with Eq. 1)21)1. where to/s can be neglected and with Ac = 0.72 ± 0.02. 

Both FigQ and FigJ2] show the accuracy of the numerical method and that the Gaussian 
approximation gives the correct results for C(t, s) and R(t, s) in d = 4 and T = Tq. Further- 
more, according to universality, the numerical values for Tc, Ar and Ac yield an amplitude 
ratio in agreement, within numerical uncertainty, with the Gaussian result TcAr/Ac = 1/2. 
No logarithmic corrections are numerically detected in the range of times investigated, as 
shown in the insets of FigQand FigGJ 

B. d = 2,T = T c 

We consider a square lattice with N = 1000 2 Ising spins and we compute numerically the 
response and correlation functions in the quench to Tc — 2.26918J, for five different values 
of s = 100, 200, 300, 400, 500. In FigCHand Figg|we have plotted the quantities 

g R (t,s) = (t-s) a+1 (t/sy 9 R(t,s) (37) 
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gc(t,s) = (t-s) b (t/s) 1 ~ e C(t,s) (38) 

versus t/s, with 9 = 0.38 and a = b = 0.115 taken from ref. [lflj . We find data collapse as 
expected from the RG results of Eqs. ()2pi3j) which yield 

gc{t,s) = Acfc{t/s), g R (t,s) = A R f R (t/s). (39) 

Using, next, the asymptotic conditions Hul^oo f R (x) = linx^oo fc{%) = 1 we can extract 
the amplitudes Ac = 0.78 ±0.02 and A R = 0.071 ±0.002. Lastly, from FigOlwe can make a 
check on the LSI prediction © that g R (t, s) ought to be constant with g R (t, s) = A R . Fig|3] 
shows that there is an evident deviation from LSI for x < 5, while the LSI behavior holds 
for x > 5. 

C. The limit fluctuation dissipation ratio X^ 

We measure X^ using Eq. (fTBj) with values for Ac, TcA R and 9 estimated from the 
numerical data. In d = 4, with 9 = 0, A c = 0.72 ± 0.02 and T C A R = 0.35 ± 0.02, we find 
-^oo = 0.49 ± 0.03 in agreement with the Gaussian result Xoo = 1/2. This supports the idea 
that not only 9, but also X^ is universal |5l. Ifil. Il6| . 

In d = 2, we find A c = 0.78 ± 0.02, T C A R = 0.161 ± 0.001 and taking 9 = 0.38 from 
Ref. [l9j] we obtain X^ = 0.33±0.01, in agreement with previous numerical results obtained 



with different methods [24J, |31| and with X^ = 0.30 ± 0.05 from the two loop e expansion 
|l8| . For convenience, the numerical values of exponents, amplitude ratio and Xoo in the 
different processes have been collected in Table HI 

D. d = 2,T<T c 

In the quench to below Tc the behavior of the data in the short time separation regime 
t-s<s allows to discriminate between the additive and the multiplicative forms of R(t, s). 
Expanding up to first order in (t — s)/s, in the former case from from Eqs. (J3J) and (jSJ) one 
ontains 



R(t,s) = R st (t-s) + s- {1+a) 



h R (l) + ti R (l)' 



(40) 
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while from the LSI form 

R L si(t,s)=A R (t-s)-<- 1+a) x° (41) 

one gets 



R LSI (t,s) =A R (t-s)~( 1+a ) 



B- 1 



(42) 



s 

Therefore, as t — s becomes small with finite s, from Eq. (f4"U|) there remains an s dependence 
due to s~ ( - 1+a ' ) /i^(l), while from Eq. (|42J1 there is no residual s dependence. 

In order to see which of the two behaviors is actually realized in the data, we have 
quenched a system of 1000 2 Ising spins to the temperature T = 1.5J ~ 0.69Tc, taking the 
wide range of s 6 [101, 1577] and focusing on the regime t — s < s. The numerical data are 
displayed in Fig0 Furthermore, in the inset we have plotted Eq. (J41J1 in the same range of 
s and t — s, with A R = 0.01 obtained by imposing R(s + 1, s) = Rlsi(s + 1, s) for s = 100, 
a = 0.27 extracted from the data for R(t, s) (see below) and 9 = X/z — 1 — a with X/z = 0.625 
|l4j |. The numerical data display an evident dependence on s down to t — s = 1, which is 
absent in those for Ri,si{t,s) (note the same vertical scale). Therefore, the LSI form of 
R(t, s) can be ruled out. 

We have also extracted R st (t — s) from the data using the following protocol. We have let 
the system to evolve in contact with the thermal reservoir at the temperature T = 1.5 J after 
preparing it in a completely ordered state, for instance all spins up. The equilibration time 
t eq for this process is finite and R s t(t — s) is obtained by measuring the response function 
for s > t eq . The data obtained in this way yield, as expected, an exponentially decaying 
contribution (continous line in FigEJ), which becomes very rapidly negligible with respect to 
the full R(t, s). Therefore, in the observed range of s and t — s, i) aging is well developed in 
the data and is due to the R ag (t, s) contribution in Eq. (j3j), while it is practically unobservable 
in the LSI and ii) the stationary contribution from the data decays exponentially, while in 
the LSI there is a power law decay. 

Next, we have extracted R ag (t, s) via the subtraction R ag (t, s) = R(t, s) — R s t(t — s) and 
we have made the comparison with the fitting formula ()25|) . which in the short time regime 
reads 



R ag (t,s)=A R s- L/z (t-s + t 



Vtft-x -J-^-H- 1 /*-« 



o A 



(43) 



s 

and predicts TTI behavior if s l l z R ag (t, s) is plotted against t — s. Indeed, this is observed 
in FigEl where we have used the exponent \j z = 0.47 J22| obtained from the numerical 
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a 


9 


T C A R 
Ac 


-^oo 


Gaussian model 


(d-2)/2 





1/2 


1/2 


Ising critical d = 4 


1.01 ±0.01 


0.00 ± 0.02 


0.49 ± 0.03 


0.49 ± 0.03 


Ising critical d = 2 


0.115 


0.38 


0.20 ±0.01 


0.33 ±0.01 


Ising T <T C d = 2 


0.27 ±0.02 










TABLE I: Exponents, amplitude ratio and Xoo for quenches to and to below Tq- The values of a 

from ref. 19J. For Xqo = in the 



and 9, in the d = 2 critical quench of the Ising model, are taken 
quench below Tq see, for instance, ref. 



data for the interface density pi(s) ~ s 1 '*. The curves for different values of s collapse on 



a master curve which is very well fitted by the power law (t — s + to) 
FigEJ • The comparison with Eq. (J43JI . then, gives the numerical value 



-0.80 



(broken line in 



a = 0.27 ±0.02 (44) 

in agreement with previous results for this exponent J3, |28[. The tiny deviations from 
the fitting curve, observed in FigEI when s is small and t — s < 2, can be attributed 
to the absence of a sharp separation between bulk and interface fluctuations in this time 
regime. This implies that Eq. (@J) is not exact for small s and, therefore, that the aging 
contribution in the response function cannot be obtained simply by subtracting R s t(t — s) 
from R(t, s). However, this procedure becomes exact for larger values of s, as demonstrated 
by the fast convergence of the numerical data for R ag {t, s) towards the behavior of Eq. (|43|) . 
Furthermore, as remarked above, the equilibrium response R st (t — s) is a very fast decreasing 
function of t — s and, when t — s > 2, the condition R s t(t, s) <C R(t, s) is fulfilled. Hence, 
even for small values of s, in the time region t — s > 2, one has always R(t, s) ~ R a g{t, s) 
and the numerical curves follow Eq. 



V. CONCLUSIONS 

We have investigated the suggestion put forward in Ref. |5|, [19j that LSI applies when 
Gaussian behavior holds, by looking at the scaling behavior of R(t, s) in the kinetic Ising 
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model with Glauber dynamics in two revealing test cases i) in the quench to Tq with d = 4, 
where deviations from Gaussian behavior are expected to disappear and ii) in the quench 
to Tq and to below Tq with d = 2 where, conversely, corrections to Gaussian behavior 
are expected to become quite sizable. Unlike Pleimling and Gambassi, who work with an 
intermediate integrated response function, we have computed directly R(t,s), in the sense 
specified in section UTTl after Eg. (|36|). producing high precision data by means of the new 
numerical algorithm of ref. J22J. In the d = 4 numerical simulation we have not found 
logarithmic corrections to the Gaussian behavior (see the insets of Fig.l and Fig. 2). 

Our results do confirm the conjecture that a Gaussian approximation is inherent to LSI. 
In the case of the d = 4 critical quench we find agreement between LSI, Gaussian behavior, 
and numerical data. Instead, in the case of the critical quench in d = 2, deviations from LSI 
are observed, which go in the same direction as the field theoretical calculations and pre- 
vious numerical results from the global integrated response functions. Similarly, important 
deviations from LSI behavior are found in the quench to below Tc- In the latter case the 
data i) are incompatible with the multiplicative form of R(t, s) predicted by LSI and ii) do 
confirm the result a = 0.27 ± 0.02 for the scaling exponent of R ag {t, s), first obtained from 
the measurements of x{t, t w ) [9J. It ought to be mentioned that the behavior of R ag (t, s) in 
the quench to below Tn of the d = 2 Ising model has already been investigated numerically 
in great detail in ref. |28(. In that paper we have produced evidence for the existence of a 
strong correction to scaling, next to the leading term behaving as in Eq. (jHJ. In the present 
work there has been no need to worry about the correction to scaling, since we have focused 
on the time sector with t — s < s and s sufficiently large, where the correcting term is 

"t::,L M also b e — U- — H. M md eon ab „ rat0ra QQ „ 

proposed a more general version of the LSI by replacing Fr(x) in Eq. JT} with 

F R (x) = A R (x- l)-( 1+a ') x e+a '~ a (45) 

where a' is new exponent. The old LSI is contained in the new one as the particular case 
corresponding to a' = a. Fitting the numerical data for the integrated response function in 
the critical quench of the d = 2 Ising model |3j|, an improvement over the old LSI has been 
obtained with a — a' = 0.187. One of the problems with the new LSI, however, is that when 
a y^ a' the numerical improvement is obtained at the expense of destroying quasi stationarity 
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in the short time regime, which is required by the separation of the time scales |l|. A detailed 
analysis of the new LSI is beyond the scope of the present work and is deferred to a future 
publication. Considerations similar to ours have been made by Hinrichsen [36j in comparing 
numerical data with the predictions of LSI J37( for the 1 + 1-dimensional contact process. 
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FIG. 1: (Color on line) R(t, s) in the d = 4 Ising model quenched to T = Tc — 6.68J. The broken 
line is the analytical solution of the Gaussian model (Eq.(|22|)) with Ar = 0.35/Tq and to = 0.1. 
Inset: plot of (t — s + to) 2 R(t,s) showing the absence of corrections to Gaussian behavior. The 
broken line indicates the constant value. 
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FIG. 2: (Color on line) C(t, s) in the d = 4 Ising model quenched to T = Tc — 6.68 J. The broken 
line is the analytical solution of the Gaussian model fEa.(|21|) with Ac = 0.72 and to = 0.1. 
Inset: plot of (x — 1 + to/s)(x + 1 + to/s)sC(t, s) showing the absence of corrections to Gaussian 
behavior. The broken line indicates the constant value. 
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FIG. 3: (Color on line) g R {t,s) denned in Eq.^ with 6 = 0.38 and a = 0.115 in the the d = 2 
Ising model quenched to T = Tq — 2.26918 J. The broken line is the prediction of LSI. 
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FIG. 4: (Color on line) gc(t,s) defined in Eq.® with 6 = 0.38 and a = 0.115 in the the d = 2 
Ising model quenched to T = Tq — 2.26918J. The broken line corresponds to the amplitude 
A c = 0.78. 
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FIG. 5: (Color on line) R(t, s) in the quench of the Ising model to T = 1.5 J below Tq. Inset: plot 
of RLsi(t, s) for the same values of s and t — s. The continous line is the plot of R s t(t — s). 
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FIG. 6: (Color on line) s 1 ' z R ag (t, s) with z = 0.47 for the d = 2 Ising model quenched to T = 1.5J 
below Tq- The broken line represents the power law behavior (t — s + £o)~ ' 80 of s 1 < z R ag (t, s) from 
Eq.@3). 
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